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As a first step toward understanding a lanscape of vacua in a tlieory of non-linear massive gravity, 
we consider a landscape of a single scalar field and study tunneling between a pair of adjacent 
vacua. We study the Hawking-Moss (HM) instanton that sits at a local maximum of the potential, 
and evaluate the dependence of the tunneling rate on the parameters of the theory. It is found 
that provided with the same physical HM Hubble parameter Hhm, depending on the values of 
parameters 03 and Q4 in the action (|2.2p . the corresponding tunneling rate can be either enhanced 
or suppressed when compared to the one in the context of General Relativity (GR). Furthermore, 
we find the constraint on the ratio of the physical Hubble parameter to the fiducial one, which 
constrains the form of potential. This result is in sharp contrast to GR where there is no bound on 
the minimum value of the potential. 



I. INTRODUCTION 



It is of great interest to promote Einstein's gravity theory to a massive one [ij . However, it has been shown that 
. Lorentz-invariant massive gravity theories would suffer from the Boulware-Deser (BD) ghost instability [H, Q 
(For a review, see e.g. Ref. Recently, a non-linear construction of a Lorentz-invariant massive gravity 

theory has been proposed in Refs. mf -lsl, where the nonlinear terms are designed so that the BD ghost is removed 
by using the constraint equations [9|, [lOj . Many researches have been made to explore its physical consequences 
since this breakthrough. 

One of the most interesting results among these explorations is that such a massive gravity theory allows 
self-accelerating solutions [ll-14|. While cosmological solutions in usual massive gravity theories are known 



to suffer from appearance of a ghost within 5 degrees of freedom of the massive graviton [15l - ll7l | , it has been 
proved that perturbations on the solution found in Ref. [l^l do not have massive degrees of freedom and evade 
the ghost instability at linear order [l^. ^ 

Though the mass of graviton may successfully explain the current acceleration, it is still to be explained 
why the cosmological constant is so small while large quantum corrections are expected. ^ In spite of many 
. attempts, this problem, so called "the cosmological constant problem" 01 HH, has not been solved yet. 
\^ ' A possible resolution to this problem is the anthropic selection of the cosmological constant in the landscape 

of vacua [20l - [2^ . In the landscape of vacua, a vacuum can be unstable with respect to tunneling to other vacua. 
The main interest of this paper is to investigate how the stability of a vacuum is determined in the non-linear 
, massive gravity theory. We investigate how the tunneling rate in general relativity (GR) is modified by the 
graviton mass terms and clarify dependence on the model parameters. Moreover, an interesting possibility will 
arise if the graviton mass depends on a value of a scalar field [1^, e.g., the tunneling field. In this case, the 
effective cosmological constant of a vacuum can be larger than the other while its potential energy is smaller 
' because of the contributions from the mass terms. It will be interesting to sec how the tunneling process 
H . proceeds in this case. 

In this paper, we consider the Hawking-Moss (HM) solution [23| for a scalar field with minimal coupling to 
gravity to understand the tunneling process in non-linear massive gravity. We set up the model and found a 
bounce solution corresponding to a HM solution. Based on this solution, we evaluate the HM action and the 
contributions from the graviton mass terms. Wc find they may either enhance or suppress the tunneling rate 
compared to that in GR. 

This paper is organized as follows. In Sec. |TT1 wc setup the Lagrangian for our model. In Sec. IIIIl wc 
formulate the equations of motion (EOM) and solve the constraint equation. In Sec. IIVI the HM solution is 
studied. To analyze the effect of the graviton mass terms in details, we consider several specific combinations 
of the parameters. 



* E-mail address: yingli@yukawa.kyoto-u.ac.jp 
t E-mail address: rsaito@yukawa.kyoto-u.ac.jp 
t E-mail address: misao@yukawa.kyoto-u.ac.jp 

^ See also the recent discussion on the ghost instability at nonlinear level [T^ . 

^ Note that the graviton mass does not receive a large quantum correction because the other parts, where the general coordinate 
invariance symmetry is respected, do not induce a correction to the mass terms. 
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FIG. 1: Illustration of the potential V{a) with two local minima of and ctt, which correspond to the false and true 
vacuum respectively. Here the atop labels its local maximum value. 



Throughout the paper, the Lorentzian metric signature is set to be (— , +, +, +). 



II. SETUP 



We study a tunneling of a minimally coupled scalar field a between two vacua in non-hnear massive gravity. 
It is described by a 4-dimensional metric 5^1/ , a fiducial metric Gab, and the Stiickelberg fields 0°. ^ The action 
is given by ^ 



where 



and 



R 



-(dar + V{a) 



C2 
C3 



- ([/C]^ - 3 [/C] [/C2] + 2 {K?\ 



— [[ICf - 6 [ICr [IC^] + 3 [IC^] + 8 [/C] [K.^] - 6 [JC^ 



(2.1) 
(2.2) 

(2.3) 



(2.4) 



(2.5) 



In this paper we assume Gab is non-dynamical as discussed below. The potential V{<7) is assumed to have two 
local minima ap and (tt, where the former corresponds to the false vacuum, and a local maximum between 
them, a = crtop (sec Fig. [T]). 



^ We use Greek letters ^,v, ... for the spacetime indices and the Latin letters ... for the space indices, while the Latin indices 
a,b,... for the internal space (Lorcntz frame) indices. Repeated indices are understood to be summed over unless otherwise 
stated. 

* We use the natural units throughout this paper. 
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Provided the action (|2.ip . the Euchdean action can be obtained by Se = iS[x'^ — > ix'^]. Then, in the 
scmiclassical hmit, the tunnehng rate per unit time per unit volume is given in terms of the Euchdean action as 

T/V = Aer^ ; B = SE[gi_i.iy,B, 4>b] - SE[gi_Ly,F,(t>F] , (2.6) 

where {5^ii/,_b, 0s} is the so-cahed bounce solution, or an instanton, a solution of the Euclidean equations of 
motion with appropriate boundary conditions, and {g^^i/.F, 4>f} is the solution staying at the false vacuum [25| . 
Among the bounce solutions, the one that gives the least action dominates the tunneling. 

Conventionally a bounce solution {g^u.Bi4>B} is explored assuming 0(4)-symmetry. This is because it was 
proved that an 0(4)-symmetric solution gives the lowest action for a wide class of scalar- field theories [l^l- It 
is therefore reasonable to assume the same even in the presence of gravity [1^ . Hence the metric would take 
the form, 

g^^dx^da;'' = N{£,fA^^ + a{£,fVl,jdx\\x\ (2.7) 
where the second term VlijAx'^dx^ is the metric on a three-sphere [K > 0), 



n = s,, + ""f ; . (2.8) 

^ 1 - KSlmX^X"^ ^ ' 

However, in massive gravity, this may not be always true if the fiducial metric Gab does not respect the symmetry. 

In this paper, to avoid possible complications due to the nature of the fiducial metric, we assume that it is 
given by the de Sitter metric with the "Hubble parameter" F, namely, 

GaMdrd(f>'' = -{d<py + b{(t>yn,jd^'d(t>', (2.9) 

where 

6(0") = F-^VKcosh{F(f)°) . (2.10) 
Then an 0(4)-symmetric solution can be obtained by setting 

0° = /(e), 0' = ^v (2.11) 

The absence of the BD ghost has been proved for a general fiducial metric in [l3] and the de Sitter fiducial 
metric has been recently investigated in [2^ [2^. In contrast to the Minkowski fiducial metric, the de Sitter 
metric is shown to lead to a wider range of cosmological solutions Q . 

At this level, the origin of the fiducial metric is unclear: It may be a non-dynamical metric of the Stiickelberg 
field space or one of the two dynamical metrics in bimctric gravity [30l - l3^ . In the later case, it may be natural to 
use the Euclidean signature also for the fiducial metric, which can be accomplished by making the replacement, 
0° — (/ — ifE), in Eq. (|2.11l) . In this case, the 0(4)-symmctric solutions would exist also in the case that 
the fiducial metric is given by the Minkowski metric and the bounce solution could be found in a similar way 
by setting 6(0°) ee ^/X0°. 

Here, we stick to the Lorentzian signature for the fiducial metric because it is non-dynamical in our case. 
Nevertheless, thanks to the assumption of the de Sitter fiducial metric, we may adopt the 0(4)-ansatz. 



III. BASIC EQUATIONS 

Here we derive basic equations. First we write down the Euclidean version of the action (j2.1|) . Inserting the 
0(4)-ansatz (|2.7p and (|2.1ip into the Euclidean version of (|2.5p . we find 



K = l-^^, ^=(l-^]5], Kl^Q, K^^Q, (3.1) 



where the dot means a derivative with respect to the radial coordinate, ' = d/d^. This gives 



[/C] = 1 



[/C^] = 1 



[IC^] = 1 - 
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Then the gravity action is reduced to 



IgE = J d^E 



-3KNa 



N 3 



ml {L2E + a^LzE + aiLiE) 



where 



C2E = 3a(a - h) 2Na - ^J-Pa - Nb 



C3E = (a - by 4Na - SyJ-pa ~Nb], 



and the action for the tunneling field is reduced to 
which coincides with the Euclidean version of the action (2.8) in Rcf. [l^ if we replace / by h{f). 



(3.2) 
(3.3) 
(3.4) 
(3.5) 

(3.6) 

(3.7a) 
(3.7b) 
(3.7c) 

(3.8) 



A. Constraint equation from the Stiickelberg field 

The introduction of the Stiickelberg field gives rise to an additional constrait equation. Variation of the action 
p.6p with respect to / gives 



(m + 7V6j) 



where 



a I \ a I \ a I \ a 



db 



= 0, 



b f= — 

df 



Ksmh{Ff). 



(3.9) 



(3.10) 



Note that the derivation of the above equation does not depends on the choice of a branch of \/ —p 
Solving Eq. ()3.9p . we obtain two branches: 

- Branch I 



Nb 



(3.11) 
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Branch II 



3 - — ) + ag f 1 - - ) f 3 - - ) + "4 f 1 - - I - . 
a I \ a I \ a I \ a, 



The solution to this equation is given by 



b = X±a , X± 



1 + 2 as + a4 ± -^/l + q;3 + q;§ — a4 



(3.12) 



(3.13) 



Hereafter, we call the choice X+ as Branch 11+ and X_ as Branch II_. In Appendix it is shown that 
Branch I is equivalent to Branch II with a different set of parameters in the case of the HM solution. Hence, 
we concentrate on the Branch II solutions hereafter. 



B. Euclidean equations of motion 

Variation of the action p.6p and p.8p with respect to N gives 



where 



Pa 



-ml \ 1 



3 2 



a3 1 



"4 1 



Inserting the solution p.l3p into the above, pg reduces to a cosmological constant, 



Pg(6/a = X±) = A± = -- 



[a-i + 014) 

Variation of the action p.8|) with respect to a gives 

1 



(1 + aa) (2 + aa + 2 - 3 04) ± 2 (l + 03 + - 



■i/2 



Na^ \ N 

Introducing the proper radial coordinate t = J Nd^, the equations are rewritten as 

3 fdaV 3K 1 fdaV . ^ 



V dr 

-12 



2 V dr 



d^o- ^ / da\ da , 



(3.14) 
(3.15) 

(3.16) 

(3.17) 

(3.18) 
(3.19) 



As explained in the previous section, the tunneling rate p.6|) is given by the Euclidean action evaluated for a 
solution of Eqs. p.l8p and p.l9p . In the following, we construct a HM solution and then evaluate the action 
for it. 



IV. HAWKING-MOSS SOLUTION 

A. Evaluation of the tunneling rate 

A HM solution can be found by setting the tunneling field to the local maximum value, o'(^) = trtop- Then 
the equation of motion p.l9p is trivially satisfied and the Euclidean Friedmann equation p.l8p reduces to 
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Setting the boundary condition as aHM(^^HM''" = '^'^1'^) = ^ind assuming A±_cff > 0, the HM solution is 
obtained as 

auyiir) = H:^1^\/K cos {H^mt) ■ (4.2) 
Here, we have introduced the Hubble parameter of the physical metric by 



Hhm^^%^ = ^^^^^4±^. (4.3) 

Inserting Eq. ()4.1|) into the Euclidian action given by Eqs. p.6|) and p.Sp . and using N^^f — d//dr, the total 
action can be expressed as 



//•Tr/2H-HM I 
d^xVO. I At a^j^j 2A±,off 

J-TT/2HnM \ 



where, for brevity, we have introduced the parameter Y± in terms of X± as 

Y± = 3(1 - X±) + 3^3(1 - X±)^ + 04(1 - X±)3 . (4.5) 

We also need the solution for the Stiickelberg field / to evaluate the action. It is given in terms of chm as in 
Eq. (13131), 

^HM = F- W cosh(F/HM) = X±aHM ■ (4.6) 



This can be solved for / as 

/hm(t) = 
where the parameter a is defined as 



"HM COs(i?HMT) ± V "hM C0s2(iJHMT) - 1 



(4.7) 



F 

aHM^X±— — . (4.8) 

As clear from the above definition, the parameter auM represents the ratio of the Hubble parameter of the 
fiducial metric and that of the physical metric. Since X± should be positive for the constraint (|3.13l) to be 
satisfied, aHM is also positive. 

Finally, we evaluate the derivative of / since the action l|4.4p depends on the Stiickelberg field only through 
dj/dr. Taking a derivative of Eq. (|4.6p with respect to r, one obtains 



d./HM X± sin(i?HMT) 



dr sinh(F/HM(T)) ' 

where sinh(F//fjv/) can be calculated from Eq. (|4.6p as 



(4.9) 



sinh(F/HM) = ± VaftM cos2(i/HMr) - 1 . (4.10) 

Though we have obtained "±" branches for d/ / dr here, the sign is irrelevant because the action (|4.4p depends 
only on its square, 

f dfnuY _ X^sin^iHiiMT) ^^^^^ 



\ dr / Q;^]y[C0s2(iJHM-r) - 1 

Here we note that wc can evaluate d/^/dr {/e = —if) in a similar way for the Minkowski fiducial metric by 
setting fofjM = \^.fE ill Eq. (|4.6|) . Taking a derivative with respect to t, one obtains 

— ^ = -X± sm{Hiiur) ■ (4.12) 

dr 

The square of this equation coincides with Eq. (|4.1ip with afjM ~ after making the replacement /e = —if- 
Hence, the analysis below includes the results for the Minkowski fiducial metric. 

We also note that the induced metric Gabdfi4>'^d,j(j)'' in the case of anM > 1 has a "singularity", where / 
becomes infinite, and becomes neither pure imaginary nor pure real. The singularity would also appear even if 
we had started from the Euclidean signature for the fiducial metric. In this case, the four sphere for the fiducial 
metric would crunch before the physical one does. For this reason, we concentrate on the case where anM < 1- 

As easily seen from the solution (j4.7p . the Stiickelberg field becomes complex in general. However, as we will 
see, the action turns out to be real for q;hm < 1- Wc also note that if we start from the Euclidean fiducial 
metric, /hm and d/nM/dr become real for ohm < 1- 
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1. Action for the case a < 1 

In the action (|3.6p . the only non-trivial term is the one that depends on d//dT. Hence, we first evaluate this 
part. Substituting the solution (|4.1ip . we find 



7r/2//HM 

-tt/2Hhm 



d/] 



HM 



dr 



-"hm 



-it/2 



cos^{Humt) sin(iJHMT) 



o„,4 rr4 
"-■"HM^HM 

-J^HM^HM 



A'5^1-22(2 + z2)|0 (z = OHM COs(iJHMT)) 



l-"HM(2 + a 



(4.13) 



As claimed, the action is real. 



2. Hawking-Moss action 



From the result of the previous subsection, we obtain the Euclidean action as 

2 



Shm = 'S'_E[aHM, CTtop] = ^TTT" 



6 v^hm; 



2 

^("hm) 



.M.i^^izzei^, ,4,4, 



where / = 2Tr^K-^^^ has been used. As previously explained, we consider only the case anM < 1 since a 

singularity would appear otherwise. In addition to the standard first term determined by the Hubble parameter 
Hum (which nevertheless contains a contribution of the mass term A±), there apprears a mass-dependent term. 
The presence of this second term may be regarded as the genuine effect of the graviton mass. We also note here 
that the function A{a) is regular at a = 0, which corresponds to the case of the Minkowski fiducial metric. 

As given by Eq. p.6p . the tunneling rate is determined by the difference between the HM action Shm given 
by Eq. (|4.14p and the false vacuum action Sp = SE{gii,^,F, 4>f)- The latter is given similarly as Eq. (|4.14p by 
replacing the Hubble parameter by 



H^^^yj^f^, ,4,5, 

and a by its false vacuum value, ap = X±F/Hp. We note that ap should not exceed unity either. Since the 
potential energy at the false vacuum is always smaller that at the maximum, V{ap) < V^(crtop), we have 

X±F 

ap > ttHM = 7j • (4.16) 

Hence the constraint ap < 1 gives a tighter constraint on the potential energy than that from Ohm < 1- 
Conversely, if the false vacuum exists without a singularity in the Stiickelberg field, a regular HM solution 
always exists. 

To see the effect of the second term, let us compare the current result with the HM action in GR with the 
same value of the Hubble parameter, B^^^) = 8Tr^{—H^^ + H^^). We have 

AB . 5 - ^ l^r.X, ( _ ^) . (4.17) 

Note that the function A{q) is positive and a monotonieally increasing function of a for < a < 1- Hence 
< A(q!hm) < A{ap). Since iJ^^^ < £fp ^, this means the sign of AB is determined by the sign of Y±. Namely, 
for y± > (< 0), we have AB < (> 0), and hence the tunneling rate is enhanced (suppressed) relative to the 
case of GR. 
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FIG. 2: The sign of Y± is shown for Branch 11+ (left panel) and Branch II- (right panel) solutions. The white and yellow 
regions corresponds to 1 + Q3 + — 04 < and X± < respectively, which should be excluded since the cosmological 
solution does not exist in these regions [13, O- The green regions correspond to Y± > 0, which implies that the HM 
tunneling rate is enhanced in comparison to the GR case for the same Hubble parameter -ffHM, while the pink regions 
correspond to Y± < 0, which implies that the HM tunneling rate is suppressed. Along the solid lines (defining the 
boundary between green and pink regions), X± = 1 and hence Y± — 0. Hence the HM solution (|4.14|l reduces to that in 
GR. The dotted lines denote X± = where the solution ceases to exist. Along the dashed lines, X± diverges and thus 
defines another boundary of the solution space. 

B. Parameter dependence of the tunneling rate 

In the above, we have seen the effect of the graviton mass on the tunnchng rate is determined by the sign 
of Y±. In this subsetion, we analyze the dependence of the tunneling rate on the parameters as and 014 by 
examining the sign of Y±. 



1. General case with mrig 7^ 



As shown in the previous subsection, the change in the tunneling rate is determined by the sign of Y±^ which 
is given in terms of X± as given by Eq. (|4.5p . Let us recapitulate it here: 



1 + 2 0:3 + Q'4 ± \J 1 + as + a\ — ai 



as + 04 

1± = 3(1 - X±) + 3a3(l - X±f + 04(1 - X±f . (4.18) 

In Fig. [21 we show the sign of Y± in the parameter space {as, a4) for Branch 11+ (left panel) and Branch II_ 
(right panel) solutions. As seen from the figure, the tunneling rate is enhanced (suppressed) for Branch 11+ 
(Branch II_) in a large region of the parameter space. It should be noted again that ap < 1 gives rise to a 
constraint on V^ap). 

Next, we focus on some limiting cases and analytically investigate the parameter dependence of the tunneling 
rate. 



In this case, from Eq. (|4.18p 



2. 03 = 



X± = — (1 + a4 ± VT^i) , 1± = (1 - ^±) (2 + a4(l - X±f) , 
q;4 



(4.19) 



and from Eq. (|3.16p . 



A± = f 

q;| 



2 - 3^4 ± 2(1 - 04)5 



(4.20) 
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Note that the parameter 04 should satisfy 0:4 < 1 for X± to be reaL In the foUowing, we evaluate the HM 
action for Branch 11+ and Branch II_ separately. 

• Branch IIj. 



In this branch, we have 



— (1 + 04 + \/l - "4) 



- — (1 + Vl - oiA(\ + Q;4 + ^/l - aA 
a-r. 



and 



A+ = i 



2 - 3a4 + 2(1 - ai)^ 



(4.21) 



(4.22) 



From the expression above, one immediately finds that 04 = can be a singular point for Branch 11+ . To 
study this case, we set as = Q!4 = e <C 1 and take the limit e — >■ 0. Then one obtains 



A+(a3=a4=0) = -^ 



2 + ^ + 0(,3) 



(4.23) 



which leads to 



lim A+(a3 = = 0) = —00 . 



(4.24) 



So as = a4 = is singular in this branch, which should be excluded in the analysis. From the condition < 1, 
we obtain the constraint on V{aF) by using Eqs. ()4.8|) and ()4.2ip : 



(27I + 37I - 1) ml + 3 (7I - 74 - 2) V2 _ (274 - l)ml + 3(74 - 2)^F' 



(l-7l)^ 



(74-1)2 



where 



74 = VT^^ (> 0) , 



(4.25) 



(4.26) 



An interesting consequence of this constraint is the absence of a (meta-)stable vacuum state below this bound. 
For example, even if there is a potential minimum below this bound, tunneling to such a state may be prohibited. 
In terms of 74, Eqs. (|4.2ip can be rc-cxprcsscd as 



74 



74 



1 



74 



(74 - ly 



(4.27) 



Thus the condition X+ > requires either 74 > 2 (04 < —3) for which we have y+ > 0, or < 74 < 1 
(1 > a4 > 0) for which 1+ < 0. Therefore the tunneling rate is enhanced for a4 < --3 and suppressed for 
< a4 < 1. The range — 3 < 04 < is excluded because a solution ceases to exist there. 

In order for the solution to be a four-sphere, the effective cosmological constant must be positive definite, 
A+.cft > 0. We find this is always satisfied, since 



A, 



,cff 



> y(crF) + A+ > SF^Xl > 0. 



(4.28) 



We again note that 74 = 1 (a4 = 0) is a singular point, which coincides with the conclusion that Branch 11+ 
does not exist for ~ ~ 0. 

• Branch II_ 

In this branch, we have 



a4 ^ 
2 



a4 — vT— 04) ~ 



74 



74 



74 + 2 
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and 



2 



2 - 3q;4 - 2(1 - q;4)5 



^74 



1 



(74 + 1) 



2"^s: 



(4.30) 



where, as before, 74 = -^/l — Q!4 and the parameter a4 should not exceed unity. We note that in contrast with 



Branch 11+ , here as 



is regular, which implies a finite A+(q;3 = = 0) ~ irn^/A. 



Similar to the case of Branch 11+ , we find a constraint on V{(jp) on Branch II_ from the condition a4 < 1 as 



-27I + 87! -l)ml + 8(7! + 74 - 2)2^2 _ -(274 + l)m2 + 8(74 + 2)2^^2 



(I-7I) 



2\2 



(74 + 1)2 



(4.31) 



It is clear that the limit 74 — >■ 1 (04 0) is regular in this case of Branch II_, which is consistent with the case 
as = a4 = studied previously. We note that there is no restriction on the range of except for the condition 
that Q!4 < 1. 

From Eqs. (j4.29p . we see that y_ is always negative. Therefore, the tunneling rate is always suppressed in 
this case. 



3. Vanishing £2 term 

In this subsection, we consider a special case when £2 vanishes but £3 and £4 are non- vanishing. To achieve 
this, we take the limit rn^ — > 0, while keeping = rriga^ and A4 = Tn^a4 finite. Substituting A3 and A^ for 
as and 014 in Eq. (|8.16|) and taking — > 0, A± is calculated as 

A± = -JX^:^^ [1 ± sgn(A3)] , (4.32) 

where sgn(A3) = IA3I/A3 denotes the sign of the parameter A3. The sign function appears due to the presence 
of the square root in Eq. p.l6|) . Similarly, from Eq. (|8.13p . X± can be obtained as 

X± = l + ^-4V [1 ± «g^(A3)] . (4.33) 
A3 + A4 

Without loss of generality, we assume A3 > since the flip of sign only exchanges the roles of Branch 11+ and 
Branch II_ . 

• Branch 11+ 

In this branch, from Eqs. (j4.82p and (|4.88|) . we have 

4:A^ 2A 

and 

m2y+ = 8^3(1 - x+)^ + ^4(1 - x+f 

= (1 - X+fX+A3 . (4.35) 

Thus m'^Y_ is always positive provided X_ > 0. From Eq. (|4.84[) . this is guaranteed if the parameter A4 is 
either in the range A4 < —8^3 or —A3 < A4. The tunneling rate is always enhanced in this case. 
Next let us consider the condition ap < 1. From Eq. (|4.34p . we find 



^ (A3 + A4)2l/(aF)-4A3 

Hence the constraint on V{<7-p) is found as 



^^'i^^^ + ^^r <1. (4.36) 



3f2(3A3 + A4)^+4Ai 
iA^TA4f • 
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Again the positivity of the effective cosmological constant, A+ oH > 0, is guaranteed because 

A+,eff > Viap) + A+ > SF^Xl > . (4.38) 



• Branch II_ 



In this branch, from Eqs. (j4.32p and (I4.33[) . wc find 

A_ = , A:_ = 1 , mlY- = . (4.39) 

This implies that the action in this case reduces to that of GR. There is no difference from GR in this branch. 
From Eq. (|4.39p . the condition ap < 1 is expressed as 

V(crF) 

which leads to the constraint on V{<jf) as 

F(ctf) > SF^ . (4.41) 



V. SUMMARY AND DISCUSSION 



In this paper, we investigated the stability of a vacuum in the landscape of vacua in a theory of non-linear 
massive gravity. For this purpose, wc derived the Hawking-Moss (HM) solution in a simple case where a 
tunneling scalar field is minimally coupled with gravity and evaluated its Euclidean action. We obtained three 
branches of the solution, which we labeled as Branch I and Branch II±, where Branch ll± corresponds to 
the self-accelerating branch in Ref . [H, [13 • We analyzed the contribution of the graviton mass terms to the 
Euclidean action, hence to the tunneling rate. We focused on Branch-IIj- HM solutions because a Branch-I HM 
solution is shown to be equivalent to a Branch-II solution with a different set of the model parameters. 

The Euclidean action of the HM instanton is found to have two distinct terms. In addition to the standard 
contribution determined by the Hubble parameter of the HM solution, we obtained a mass-dependent non- 
standard term. To study the effect of this latter non-standard contributions, we compared the tunneling rate 
to that in GR for the same value of the Hubble parameter. 

For Branch IIj-, the enhancement or suppression of the tunneling rate relative to GR is determined by the 
sign of a quantity, denoted by Y±, and it depends on the model parameters and a^. We found that in a wide 
area of the parameter space the tunneling rate is enhanced for Branch 11+ and suppressed for Branch II„ . 

Wc also found that the solution can exist only if the Hubble parameter of the ratio of the physical metric to 
the fiducial one is greater than X±. As a result, the form of the potential is constrained. This seems to imply 
that tunneling from or to a vacuum whose energy density is less than a critical value is prohibited. This result 
is in sharp contrast to GR for which there is no bound on the minimum value of the potential. 

As tunneling mediated by a HM instanton is considered to be equivalent to stochastic process of going over 
the barrier due to large vacuum fluctuations intrinsic to light scalar fields in de Sitter space [36j , it is interesting 
to understand the role of the graviton mass terms in the context of stochastic dynamics. 

We found that the tunneling rate can be enhanced or suppressed depending on the model parameters for a 
simple case of a tunneling field minimally couples with gravity. Though we discussed this simple case as a first 
step, it will be interesting to see how the tunneling process proceeds when the model parameters such as the 
graviton mass depend on the tunneling field. In this case, the effective cosmological constant of a vacuum can 
be larger than the other while its potential energy is smaller because of the tunneling field dependence of the 
graviton mass. 



Appendix A: HM solution for Branch I 

In the text, we focused on HM solutions based on Branch II constraint p.9|) . Here, we show that the HM 
solution in Branch I is equivalent to the one in Branch II with a different set of the model parameters. 
Inserting Eq. (P?TU| into Eq. (|XTT|) . we obtain 



F^b^ + a'^ = K . 



(Al) 
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Hence 



Inserting this into Eq. p.l4p and setting a = crtop, we obtain an algebraic equation, 
where 



Pg = -771^ (1-61) 3 {2 -9) + as (1 -0){A-d)+ (1 - 6] 
Hence, d — hja becomes constant. The HM solution in Branch I is obtained as 

"hmW = {eFy^yfK cos {OFt) . 



F^ 



b 
a 



(A2) 



(A3) 



(A4) 



(A5) 

Comparing this with Eqs. (|3.13p and (|4.2p . we see that the form of the solutions are the same except that X± 
and TJhm are replaced by 9 and 9F, respectively. We note that the quantity corresponding to a = X±F/H-tiM 
in Branch I is unity. 

Thus, we conclude that the analysis of the HM solutions for Branch II can be directly applied to those for 
Branch I if we make the replacement X± — >■ 6 and i?HM — > OF. Note that 9 depends on F, rUg, and V{atop), 
as well as and 0:4. Hence, the dependence on the model parameters becomes more complicated in Branch I 
than in Branch II. 



Appendix B: List of symbols 



Throughout the paper, there appear many parameters which are denoted by various symbols. For convenience, 
we list some of them who play important roles in the table HI 





TABLE I: List of s 


ymbols 


Symbols 


Definition 


First appearance 


x± 


SeeEq. (|3.13|) 


(13.131 


A± 


SeeEq. (|3.l6|) 


(|3.16| 


A±,eff 


V^(o-top) + A± 




y± 


See Eq. gS)) 




Hum 


VA±,eff/3 


631) 


Qhm 


FX±/Hhm 


dm 


Hf 


^{V{aF) + A±)/3 


l|4.15|l 




FX±/Hf 


Below (14.151) 
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